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q\ • Abstract 

We present a recursive method to calculate a large q expansion of the 2d q- 

^ \ states Potts model free energies based on the Fortuin-Kasteleyn representation of 

the model. With this procedure, we compute directly the ordered phase partition 

function up to order 10 in 1/y/q. The energy cumulants at the transition can 

be obtained with suitable resummation and come out large for q < 15. As a 

consequence, expansions of the free energies around the transition temperature 

are useless for not large enough values of q. In particular the pure phase specific 
i — i 

heats are predicted to be much larger, at q < 10, than the values extracted from 

■ 

Q\ • current finite size scaling analysis of extrema, whereas they agree very well with 



recent values extracted at the transition point. 



O-t, Resume 

<D . 

Une methode recursive pour calculer un developpement a grand q du mo dele 
de Potts bi-dimensionnel a q etats est presentee, sur la base de la representation 

de Fortuin-Kasteleyn. Avec cette procedure la fonction de partition dans la phase 
' x 

ordonnee est calculee directement a l'ordre 10 en l/y/q. Pour q < 15, les cumu- 
lants de l'energie sont trop importants pour rendre utilisable le developpement 
de l'energie libre au voisinage du point de transition. En particulier les chaleurs 
specifiques predites pour les phases pures sont beaucoup plus grandes que les 
valeurs extraites des analyses de taille finie d'extrema pour q < 10, alors qu'elles 
sont en tres bon accord avec celles recemment obtenues au point de transition. 
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1. Introduction 



The two-dimensional (/-states Potts model is a very useful framework for prob- 
ing numerical algorithms and methods to analyze transitions. The advantage of 
this model is twofold. First many of its properties, such as the location of its 
transition temperature Z^ -1 and the internal energy densities in the first order 
transition regime, are known [1] and put strong constraints on the analysis of 
any numerical results. Next, its simplicity makes it easy to implement and the 
existence of the free parameter q (the number of states) allows one to vary the 
properties of the model, especially the magnitude of the correlation length relative 
to accessible lattice sizes. 

Two recent analytical results are of interest to our purpose. It has been 
shown [2] that close to /?t, the partition function Z of the Potts model, in a box 
of volume V = L 2 with periodic boundary conditions, is equal to the sum of 
the 'partition functions' Zi = exp(VFi) of the q ordered and of the disordered 
pure phases, up to a correction that falls off exponentially faster with L. The i th 
phase free energy Fi is V independent and differentiable many times with respect 
to the inverse temperature (3 at (3%. As a consequence of a recent calculation of 
the disordered correlation length at /?t[3], the interface tensions also have been 
analytically predicted [4] . 

On the numerical side, various simulations have recently been performed for 
different values of q between 7 and 20, with the main purpose of accumulating more 
experience on the identification, by numerical means, of the nature and properties 
of a phase transition. Although the general overview acquired seems quite con- 
sistent, there remains some unsatisfactory issues such as slight inconsistencies in 
finite size scaling analysis of the energy cumulants close to the transition tempera- 
ture /3 t _1 , and some discrepancies between exact results and numerical simulations 
for the interface tension. A proper understanding of these differences is important 
since in other cases of physical interest such as the 3D q = 3 Potts model or QCD 
at the deconfinement transition, no good analytic solution exists and one has to 
resort to similar numerical calculations to determine these quantities. 

The problem at hand can be stated as that of disentangling finite size effects 
coming from adding up the "asymptotic" partition function Zi, and those asso- 
ciated with truly non-asymptotic contributions coming from interfaces between 
coexisting phases. One way towards the solution of this problem is to learn more 
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about the pure phase free energies -Fj(/3) in the vicinity of (3 t and as a function of 

This is our motivation for starting a large q expansion of F t for f3 near (3 t 
based on the Fortuin-Kasteleyn representation [5]. Such an expansion is much 
in the spirit of the pioneering work by Ginsparg, Goldschmidt and Zuber [6], 
who pointed out that in d-dimension, at any finite order in z = (/ _1 / d , only a 
finite number of terms contribute in the Z 9 character expansion of the partition 
function. Low temperature expansions also exist for the same model [7], but the 
above method seems more adapted to our goal of studying the model close to the 
transition temperature. A short account of our work has been already published 
with the free energy expanded in z = 1/^ up to order 9 [8] or 10 [9]. The 
present paper provides a detailed description of the expansion and an analysis 
of the behavior in q of the six first energy cumulants. More phenomenological 
consequences on finite size analysis of numerical data, sketched in [8,9], are fully 
developed in a separate publication [10]. 

The organization of the paper is as follows. In section 2 we recall some useful 
properties of the 2-d (/-states Potts model and its representation as a model of 
bonds. Section 3 is devoted to the detailed computation of the partition function. 
The resulting free energy expansion is given in section 4. Resummation techniques 
based on Pade approximants are applied in section 5 to the energy cumulants at (3 t . 
Strong evidence is given that they increase very fast as q is lowered toward 4, the 
more so their order increases. Quantitative predictions are given down to q = 6, 
some of them being in clear disagreement with values extracted from numerical 
simulations not analyzed at (3f. Concluding remarks are presented in section 6. 

2. The model and its Fortuin-Kasteleyn representation. 

We consider the two-dimensional (/-state Potts model [11] on a square lattice 
with V = L 2 sites, defined by the partition function 

Z=Y,e X p{-0H), H = -J2^: (1) 

Wi} (ij) 

where i and j denote the lattice sites, (ij) the pairs of nearest neighbors and 
(7j = 1,2, ■••,(/. The symmetry group of the Hamiltonian is the permutation 
group of q elements. 

There is abundant literature on this model [1], and many of its properties are 
known exactly. In particular, it possesses a phase transition which, for q > 4, is 



2 



first order [12] and lies at a coupling 



A = log(Vg+l). (2) 

At (3 = Pf, Q + 1 phases can coexist. One is the disordered phase, the q other 
ones are degenerate ordered phases. The internal energy densities at (3 = /?t in the 
disordered and ordered phases are E d and E Q respectively, with 

E d + E Q = -2(1 + ^) (3) 

and a latent heat 

C = E d -E = 2(1 + tanh - JJ [tanh (n6)] 2 (4) 

y q 1 n>0 

where 2 cosh 9 = ^fq. Duality relates the ordered and disordered free energies 

F d (p) = F (p)-2\n((^ -l)/^q) , (e? - l)(e" - 1) = q. (5) 

A large q expansion can be obtained through the Fortuin-Kasteleyn [5] rep- 
resentation of the Potts model partition function 

Z = J2(^-l) b{X) q c{X) (6) 
x 

where X is any configuration of bonds on a square lattice, b(X) is the number 
of bonds in configuration A, and c(A) its number of clusters of sites : two sites 
bound to each other belong to the same cluster ( an isolated site is a cluster). 

The completely ordered configuration, A Q , has all possible bonds and, hence, 
corresponds to c(A Q ) = 1 and b(X Q ) = 2V. So the partition function of the ordered 
phase, Z Q can be reorganized as an expansion in q~^ about A Q : 

Z = q^-lf V KAV)C^)- l q k -K (7) 

fc>0,Z>0 

where N£ j(V) is the number of configurations in a volume V with / bonds removed 
and comprising k+ 1 clusters. We made explicit in Eq. (7) the factor (e 13 — l)/y/q 
which is 1 at the transition and thus provides an expansion of Z Q in q~ x l 2 near 
the transition. The enumeration of all the N% t (V) such that (/ — 2k) < M yields 
this expansion to order M in g -1 / 2 . 
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Similarly, the completely disordered configuration, X d , corresponds to c = V 
and 6 = 0. Hence we reorganize the expansion of the partition function of the 
disordered phase about this configuration at temperature as 

z * = <i v E N U y )(^) l ^-\ (8) 

j>o,z>o 

where N^ L (V) is the number of configurations in a volume V with / bonds and 

V — j clusters. Since the suppression of j clusters requires that at least j bonds 
are restored, one has k = I — j > 0. Hence we rewrite Zd as 

Z* = q V E Ntl(V)C-^) l Q k -K (9) 

l>0,k>0 

where t (V) is the number of configurations in a volume V with / bonds and 

V — 1 + k clusters. Duality as given in Eq. (5) implies N% t (V) = iV£j(V), a trivial 
geometrical property exemplified later on. Thus, in the following we restrict our 
calculation to the ordered phase contributions, except for comments on duality. 

With a large enough volume V (to eliminate all finite size effects) and periodic 
boundary conditions (to eliminate edge effects), N(V) are polynomials in V to 
any given finite order M, and all configurations retained correspond to disordered 
(ordered) islands in a bulk ordered (disordered) phase. 

3. Evaluation of the expansion 

This simple expansion which computes directly Z Q can be made recursive. 
To this end we explicitly remove successively up to 8 bonds in sub-section 3.1 
where the main properties are exhibited in order to make the general construction 
in sub-sections 3.2-3.6. To simplify, the number k of clusters added to the bulk 
ordered cluster will be called in the following the number of cluster. 

3.1. Explicit first steps 

Starting from the ordered configuration, we have 2 V ways to remove one bond, 
the number of clusters remaining 0, so that Nq i(V) = 2V. Continuing to remove 
bonds, we obtain N^ 2 (V) = 2V(2V - l)/2 and N§ t3 (V) = 2V(2V - 1)(2V - 2)/6. 
To proceed, let us introduce C(n, I) as the number of ways to remove / identical 
bonds among 2V — n ones 

1 Z_1 

C(n,l) = -J[(2V-n-i), (10) 

' i=0 
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Thus 



C(0,1), 
C(0,2), 



C(0,3). 




A new situation occurs with 4 bonds re- 



moved. We have C(0, 4) ways to remove them, 
but some of them give rise to an isolated 1- 
site cluster as shown in Fig. 1, when the 4 



Fig. 1 

Four removed bonds (crosses) 
making a 1-site cluster. 



removed bonds have one site in common. The number of ways to make this figure 
is V, the number of possible positions for the one site cluster. Thus we have with 
4 removed bonds 



where the number of configurations which do not correspond to the specified num- 
ber of clusters in Nq A is subtracted. 

With 5 removed bonds, nothing new happens and one has 



where the factor C(4, 1) takes into account the number of ways to remove one bond 
when 4 bonds have already been removed to make the configuration of Fig. 1. 
With 6 removed bonds, the new type of 

configuration shown in Fig. 2 can be made with 

2V possibilities (V translations and 2 orienta- >< >< 

tions), and its corresponds to one cluster (the ^ ^ 

2 sites surrounded by the removed bonds are 
in the same cluster). Thus we have 



Here again the multiplicity C(0,6) is corrected for all the 1 cluster contributions. 



K A (V) =C(0,4)-V, 



N? iB {V) =VC(4,1), 
A^OO =C(0,5)-yC(4,l), 



N? fi (V) =FC(4,2) + 2y, 
^ °, 6 (^) =C(0,6)-VC(4,2)-2V. 



Fig. 2 

Six removed bonds (crosses) mak- 
ing a 2- site cluster. 
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Next with 7 bonds removed, the higher 
contribution comes from the configuration of 
Fig. 3 with 2 clusters and with 2V possibili- 
ties, and thus N$ 7 = 2V. The contributions 
to N° 7 are obtained from the configurations 
contributing to iV° 6 with one more bond re- 
moved. First we have V C(4,3) configurations 
made of Fig. 1 with 3 more bonds removed, 
out of which we have to subtract the 4V con- 
figurations with 2 clusters as in Fig. 3 (the 3 
removed bonds can surround the 4 neighboring sites of the 1-site cluster of Fig. land 
these configurations do not give new contribution to N$ 7 for which the correct 
counting is already made). Next we have to remove one bond to the 2V config- 
urations of Fig. 2. If we remove the "internal" bond connecting the two sites of 
the cluster, we will get 2 clusters. So we only remove bonds of the bulk ordered 
phase, that is only 2V — 7 bonds can be removed. Thus we have 



Fig. 3 

Seven removed bonds (crosses) 
making two 1-site clusters. 



No >7 (V) =2V, 

N°j(V) =VC(4, 3) - 4V + 2VC(7, 1), 
N^ 7 (V) =C(0, 7) - N? t7 (V) - N2 t7 {V). 



As a last simple example, let us consider 
8 removed bonds and discuss the new contri- Fig. 4 

butions beside the contributions deduced from S'iies forbidden (circles) for a 

N£ 7 with the appropriate changes in the C second 1-site cluster. 

factors. New configurations with 2 clusters can be made by putting the Fig. 1 
twice on the lattice. Once the first figure is put on the lattice (V possibilities), the 
second one cannot be put at the same site nor at the 4 neighboring sites as shown 
by open dot in Fig. 4 because a bond ending there has already been removed. 
Thus the number of such configurations is V(V — 5)/2 with the 1/2 factor for 
symmetry. Concerning the contribution N° 8 , we first correct the contribution 
coming from Fig. 1 plus 4 removed bonds forming themselves a 1-site cluster, that 
is a correction V(V — 5) ( no 1/2 symmetry factor). Then we add the contribution 
coming from the configurations shown in Fig. 5 where the cluster extends either 
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over 3 sites at left ( 6V possibilities) or over 4 sites at right (V possibilities). Thus 
N^(V) =2VC(7,l) + lv(V-5), 

N° S (V) =VC(4, 4) - 4VC(7, 1) - V(V - 5) + 2VC(7, 2) + 6V + V, 

N^ 8 (V) =C(0,8)-AT° 8 -iV 2 V 

In N° 8 we keep the two contributions of Fig. 5 separate because at the next step, 
as we explained for N° 8 from Fig. 2 with one more removed bond, it is easiest to 
restrict the removed bond to belong to the bulk ordered phase, that is 2V — 10 
and 2V — 12 possibilities respectively. 



Fig. 5 

One cluster with 8 removed bonds, extending 
over 3 sites at left and 4 sites at right. 

These examples illustrate the way how a recursive construction of the expan- 
sion up to a given order M can be made. This construction involves 5 steps: 

i) Construct the dominant configurations made of one connected set of 1-site 
clusters (two clusters with neighboring sites belong to the same connected set 
). We note (k,l) the set of configurations with k connected 1-site clusters 
and / removed bonds and they are constructed for all the k and / such that 
I — 2k < M. Thus Fig. 1 and Fig. 3 belong respectively to (1,4) and (2,7). 

ii) Construct the sub-dominant configurations made of one connected set of 
clusters involving at least one cluster with more than one site ( cf. Figs. 2 
and 5 ). They can be obtained by restoring one or several internal bonds of 
dominant (parent) configurations and we note (k, I; k p , l p ) the configuration 
set with k clusters and I < M — 2k removed bonds inside a parent with k p 
clusters and l p removed bonds. Thus the configuration of Fig. 2 belongs to 
(1,6;2,7), those in Fig. 5 left and right to (1,8;3,10) and (1,8;4,12) respectively. 
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iii) Construct the product configurations which are disconnected sets of clus- 
ters. Once we have all the connected sets, we compute the the number of 
ways to put together on the lattice connected sets of clusters with Zj and hi 
such that Yl h ~ 2 Yl ki < M. For example the contribution V(V — 5)/2 to 
A^ 2 ,8 found above is the product (1,4) * (1,4). 

iv) Construct the correcting contributions which are obtained from the con- 
figurations with k and / of steps i) to iii) plus n removed bonds (1 < n < 
M + 2k — I) and have a number of clusters higher than k. These correcting 
contributions will allow to get from the contributions of configurations ob- 
tained in steps i) to ii) weighted by C(l,n) the correct contribution to N% l+n 
by suppressing the contributions from higher number of clusters (as the cor- 
rection AV in the contribution to N° 7 of Fig. 1 plus 3 removed bonds). 

v) Collect the results of steps i) to iv) to get all the N k /s relevant to order 
M. In the steps i) to iv) the number of removed bonds is minimal in the sense 
that they are all necessary to make the clusters. Thus this step collects the 
preceding contributions with appropriate signs (for correction) and C factors. 
These different steps are explained in details in the following sub-sections. Of 

course any configuration drawn on the lattice has V copies obtained by the lattice 
translations and two configurations will be said distinct if they are unequivalent 
upon lattice translation. 

3.2. The dominant configurations 

If we consider all the connected sets of clusters corresponding to k clusters 
with / removed bonds and extending on n given (connected) sites, the exponent 
/ — 2k of q-1/2 in Z Q is minimum for k = n. Hence the name of dominant for a 
(k,l) set of configuration extending on k sites and thus made of k 1-site clusters. 
All the dominant (/^-configurations as defined by the data of k connected sites 
can be obtained recursively from the (k — l)-ones 

i) connecting one additional site in all the possible ways to all the (k — 1) con- 
figurations, keeping only the distinct (upon translation) configurations. 

ii) eliminating those configurations which happen to have more than k clusters ( 
extra sites can be isolated from the bulk ordered ones and these configurations 
can be identified by a suitable cluster finding algorithm [13] ). 

In this construction, a dominant (k, /)-configuration contributing at order 
m = I — 2k can be obtained from one or several of the following dominant k — 1- 
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configurations 

i) (k — 1, 1 — 3) contributing to order m — 1 

ii) (fc — — 2) contributing to order m 

iii) (fc — 1, / — 1) contributing to order m + 1 

An example of case i) is given by Fig. 3 obtained from Fig. 1, while examples of 
case ii) and iii) are given in Fig. 6 left and right respectively. 




Fig. 6 

The same (6,17) dominant contribution obtained by connecting a 6 th 1-site cluster (open 
point) to a (5,15) (left, case ii) or to a (5,16) (right, case iii). 

Note that in the latter case, the order in g -1 / 2 decreases in the step k — 1— > k. 
Hopefully, the example of Fig. 6 is in fact generic, and one can convince oneself 
that a (k, ^-configuration obtained via iii) can always be obtained via either ii) or 
i). This means that, once the maximum order M of the computation is given, it 
is never necessary to keep contributions to a higher order in the iterative process. 

Once all the distinct dominant configurations with k 1-site clusters are ob- 
tained, they can easily be classified according to the number / of removed bonds. 
Thus a (k, I) set is given by a set of 2k data, the x and y positions of the k 1-site 
clusters. Their number is given in Table 1 and represents the corresponding con- 
tribution to N^t/V. At order M = 10, 59 (k,l) set are contributing, the largest 
size to be considered being a (25,50) configuration, a 5 x 5 square of 1-site clusters. 

3.3. The sub-dominant configurations 

Sub-dominant configurations correspond to connected set of clusters involving 
at least one cluster extending over more than one site. They can be obtained from 
the dominant configurations by restoring some of their 4k — I internal links. 
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Let us start from a (/c p , l p ) dominant set, the parent of a family of sub- 
dominant (fc, /; k p , l p ) configurations. If one interior bond is restored, the two sites 
it links now belong to the same (2-sites) cluster. Then k = k p — 1 and I = l p — 1 
and the exponent is one unit less that the parent one ( m = I — 2k = m p — 1 ). 
Hence the name subdominant. For example Fig. 2 can be obtained from Fig. 3 by 
restoring the interior bond (only one possibility), leading to a 2-site cluster. 

These sub-dominant configurations are obtained from the dominant (k p J p ) 
configurations by restoring the interior bonds one after the other in all possible 
distinct ways, looking at the k and / values of the generated configurations. In 
fact, it is not necessary to keep memory of the restored bonds and for each (k p , l p ) 
configuration it is enough to count the number of distinct (k, I; k p , l p ) configura- 
tions obtained. Thus a (fc, /; k p , l p ) set is given by a set of 2k p + 1 data, the x and 
y positions of the parent k p 1-site clusters and their multiplicity. 

Restoring systematically bonds among the Ak p — l p interior ones can be time 
consuming for large clusters. However restoring up to 5 bonds can be implemented 
easily because there is a classification according to the number of 4 1-site cluster 
making a plaquette, with the property that all the parent configuration give rise 
to the same multiplicity. 

For the order M = 10 considered here, we had to construct 146 sets (k, /; k p , l p ) 
from the 59 dominant (k, I). 

3.4. The product configurations 

Let us now consider disconnected configurations. They are made of com- 
ponents which are either dominant or subdominant configurations. In fact for 
components which are sub-dominant configuration, the corresponding result is the 
same as given by the 'parent' component, and thus we can only consider product 
of dominant contributions. Here "disconnected" implies that two requirements 
must be fulfilled 

i) no overlap between links of any two of the connected component; 

ii) no site of the surrounding bulk ordered phase become a cluster. 
According to these requirements, the values of k and / of a product are the 

sums of the k^s and /j's of the factors (its components). The method used is 
first to count all the possible ways the components can be put on the lattice with 
requirement i) satisfied, then subtract the number of configurations which do not 
satisfy k = J2^i- 
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For the product of two configurations given by the sets F\ and F 2 of occupied 
sites, once the first configuration Fi is put on the lattice with V possibilities, we 
have to search all the possibilities to put the second one. For that purpose let us 
define as F12 the set of sites around F\ where a given site of F 2 taken as origin 
cannot be put according to requirement i). If V12 is the number of sites in Fi 2 , 
then we have V(V — U12) configurations for the product F\ * F 2 (divide by 2 for 
identical configurations). An example is the product (1,4)*(1,4) as already been 
given in subsection 3.1. 

This process can be generalized to product of more figures. For the product 
of Fi* F 2 * F 3 let us first define the protected areas F 12: F 13 and F23 with F t j 
the set of Vij sites around Fi which cannot be occupied by an origin site of Fj 
according to requirement i). Once Fi is on the lattice (V possibilities) let us first 
put F 2 far enough from Fi such that F 3 can be put every where except on F 13 or 
F23. This will be the case when F 2 is every where except in a region -Fi 2 ,23 for 
which there will be overlapp of sites of F13 and F23. If 1*12,23 is the number of sites 
of Fi2, 23, we get a first number of possible configurations which is 

n a = V(V - Ul2,23)(^ - v 13 - v 23 ). 

Next we consider all the possible positions of F 2 inside ^12,23 but outside Fi 2 
(v 13,23 — v 12 possible positions). For each such position there is a protecting region 
F( 12 )3 with i?(i2)3 sites forbidden for F 3 according to rule i) (17(12)3 < v 13 + ^23) 
and we have to add to n a 

■^13, 23— V12 

U b = V J2 ( V - V (12)3)- 

i=l 

This method can be extended to products of more factor and implemented on 
computer. As an example the resulting polynomials in V for the powers of (1,4) 
needed at order M = 10 are 

(lA) 2 =lv(V-5), 

(1,4) 3 =-V(V 2 - 15V + 62), 

\ (11) 
(1, 4) 4 = — V(V 3 - 30V 2 + 32W - 1254), 

(1, 4) 5 = — V(V 4 - 50V 3 + 995V 2 - 9370V + 35424). 
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Finally we have to fulfill requirement ii), and thus subtract to the above 
configurations those which do not correspond to the specified number of clusters, 
Y^ki. Two examples are given in Fig. 7, with at left a peculiar configuration of 
(1,4) 4 contributing to a 5 clusters contribution and at right a product (3, 10) 2 
contributing to 7 clusters. 



Fig. 7 

Product of configurations giving rise to extra cluster number. 
These corrections can be obtained by trying to cast dominant or subdominant 
configurations with I = J2k into the corresponding products (see next subsection). 

3.5. Correcting configurations 

In the three preceding sub-sections, it was assumed that the number of re- 
moved bonds was minimum to make a definite configuration of clusters. In other 
words there is no "free" links, that is removed bonds not belonging to any cluster. 
Thus the number n^,l) °f configurations (k, I) gives directly its contribution to 

N Kl 

[contrib. (k, I) to N k ,i] = n (jM) . 

(For simplicity we consider only dominant configurations). We now wish to get its 
contribution to Nk,i+i, 1 < i < M + 2k — /. This contributions is n^ k ^C(l,i) up 
to corrections from the configurations corresponding to a higher value of cluster 
number. These corrections are twofold. First these i removed bonds can make 
a cluster and the corresponding corrections can be computed as in the preceding 
subsection. Or the configuration (/c, /) with i removed bonds can make connected 
sets of clusters which can be considered as inside dominant configurations. Thus 
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in Fig. 8 a (2,7) configurations with 5 more removed bonds can make connected 
sets of clusters inside a (4,12) at right or inside a (4,13) at left. The distinction be- 
tween the different possible parent dominant configurations is necessary for further 
removing of bonds. 




Fig. 8 

(2,7) configuration plus 5 removed bonds making connected 
sets of clusters inside a (4,12) at right or a (4,13) at left. 
If n((k,i)+i-.kP,iP) is the number of these r configurations , then 

[contrib. (k,l) to N k j +i ] = n (M) C(Z,z) - n^+^P^, 

with eventual summation on r. This process has to be continued and let us suppose 
that the configuration ((k, I) + i; k%, 1%) plus j removed bonds can make higher con- 
figurations (((&, /) + z; k%, 1%) + j; fcf, l p s ). Then there is also higher configurations 
obtained directly from (k, I) plus i + j removed bonds inside the same parents. 
The proper correction which take correctly the % + j bonds symmetrization is the 
direct one. This means that we have to correct the correction itself on the form 

[contrib. (k,l) to N k ^ +i+j ] = n {kyl) C(l,i + j) 

- [{n{{k,l)+i;k p r ,l%) C ( l ri j) - n (((k,l)+i;k?,iv)+j;kv,l p s ))] 

-W((fc,I)+ i+ j;fcJ 1 i;). 

with eventual summation on r and s. 

Clearly this computation apply to dominant or sub-dominant contribution as 
well as to product of connected sets of clusters and for all the configurations from 
subsections 3.2 to 3.4 characterized by k and /, we have in this step to compute 
recursively all the connected sets of clusters obtained with / / more removed links 
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(1 < If < M + 2k — I). The recursive process means that we start from all the 
already obtained configurations with (k, I) plus l'f removed bonds (0 < l'f < If) 
and add If — l'f removed bonds. 

Such configurations will be in the list of dominant or sub-dominant contri- 
butions. So as in section 3.3, we have to search all the possibilities to found a 
starting configuration plus the removed links inside the dominant configurations 
in such a way that all the exterior links match. In the example of a (2,7) con- 
figuration and 5 removed bonds, we know from the sub-dominant configuration 
analysis that a connected set of clusters with 12 removed bonds and more than 2 
clusters can be inside a (4,12) or a (4,13). For the (4,12) there is 2 possibilities 
to found the (2,7) configuration of Fig. 9. For the (4,13) configurations, we have 
to search the possibilities to put Fig. 9 inside the (4,13) with one restored interior 
link and found that 3 (4,13) configurations give rise to 2 possibilities and 8 other 
ones to 1 possibility. 

3.6. Collecting the contributions to N k ,i. 

In computing the partition function from Eq. (7), we need the value of Nkj 
for each values of k and I. The contributions to N k ,i comes from the dominant, 
sub-dominant and product configurations obtained in subsections 3.2-3.4 and cor- 
responding to a number k c of clusters and a number l c of removed bonds such 
that k c = k , l c < /, up to corrections as explained in subsection 3.5. The whole 
procedure is better understood on a particular example and let us consider the 
computation of the N 2j ±2 coefficient. Its contributions comes from dominant, sub- 
dominant and product configurations with k = 2 and / < 12 which can be 
. (2,7) dominant configurations, 

. 2,9), (2,10), (2,11) and (2,12) sub-dominant configurations, 
. (1,4)*(1,4), (1,4)*(1,6), (1,4)*(1,8) and (1,6)*(1,6) product configurations. 
For each of them, we have to correct their naive weight (number of configurations 
times the appropriate C factor of Eq. Eq. (10) for the extra links) for the config- 
urations corresponding to a higher number of clusters. The various contributions 
building up the final value of N 2y i2/V are given in the first column of Table 2, 
along with its corresponding set of cluster in column 2 (the last column labels the 
lines). Note that when a C factor is present, its second argument gives the number 
of free links, that is the number of removed bonds not belonging to a cluster. 

The first contribution (line 1) comes from the dominant (2,7) configurations 
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(two such configurations weighted by C(7, 5), the number of ways to removed 5 
bonds among 2V — 7 ones). Its corresponding corrections (contributions corre- 
sponding to a number of cluster higher than 2) are given in the lines 2 to 7. First 
in line 2, 4 removed bonds can isolate a 1-site cluster as in Fig. 1. For each (2,7) 
configuration there is (V — 8) such possibilities ( the product (2, 7) * (1, 4) ) with 
1 link remaining free among 2V — 11. Lines 3-6 concern clusters obtained with 
(2,7) and 3 to 5 removed bonds. We made a distinction between different possible 
parents in lines 5 and 6 because they imply different C factor for their correspond- 
ing contributions at AT 2 ,i3. Furthermore these (2,7)+5 corrections represent the 
whole contribution to higher cluster numbers obtained with (2,7) and 5 extra re- 
moved bonds such that all the 5 extra links belong to the clusters. However in the 
line 3 such corrections are included and should not have been subtracted. Line 7 
provides this "correction to the correction" of line 3, and thus with a positive sign. 

Lines 8 and 11 to 15 concerns contributions of sub-dominant contributions. 
They are characterizes by their parent distribution in order to determine the num- 
ber of frozen links which give the arguments of the C coefficients. Only line 8 
needs subtraction in line 9 and 10, either 2 or 3 extra links giving rise to k > 2 
configurations. Let us note that there is two ways to have a 2 clusters inside 
a dominant (5,15), one in line 14 when restoring 3 bonds, one in line 13 when 
restoring the four bonds of a loop. 

We finally start with The possible products of dominant or sub-dominant k = 
1 configurations ( lines 16,23,25-27). First in the product (1,4)*(1,4) we consider 
the two (1,4) as different in such a way that all the corresponding contributions 
(lines 16-22) have an explicit 1/2 factor. The first correction concern the 4 extra 
removed bonds making a (1,4) cluster. This product (1,4) 3 have no extra symmetry 
factor, the third (1,4) factor being formed by the extra links distinctly from the 
two other (identical) factors. The connected cluster corrections are in line 18-21, 
with again a "correction to the correction" with a positive sign in line 21. We have 
also to make a disconnected cluster correction in line 22, when 3 of the extra links 
are making a 2 clusters with each of the product configurations. For the other 
products only one of them need to be corrected. 

All these computations have been made automatic on workstation. 

4. Results 

We have computed the series for Z Q up to order M = 10 in <? -1 / 2 , which 
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involves contributions up to Nq 0j 2 5 (V), that is up to 60 removed bonds and 25 
clusters. This series exponentiates in V to give F Q to order M in the form 




+ 0(g- n / 2 ), (12) 



where the non-zero coefficients A^ m are given in Table 3 up to m = 10. 

One may notice, at the top of the columns, the appearance of stable sequences 
with decreasing /. For k even, the m/2 higher / values are in the sequence 1, 6, 22, 
68, 187, ... that is as long as the contributions are coming from removing a corner 
(case ii) when starting from a square of size m/2 m/2. For m odd, the series is 2, 
8, 30, 88, ... for the (m-l)/2 higher / values for the same reason when starting from 
a rectangle (m+l)/2 (m-l)/2. In contrast the bottom of each column is totally 
alternating. 

Our series truncated at / = 20 can be compared to the result of Ref. [7] . Up 
to order 9 in g -1 / 2 , the two series are in agreement. At order 10, some coefficients, 
missing in [7], cannot be compared, but we disagree with the coefficient of r^uir 
of Ref. [7] (we find 3822 instead of 3818) and of r 2 u 16 (6269 instead of 6265). 

Our result is obtained for the ordered phase. Then the disordered free energy 
can be obtained using the duality relation given in Eq. (5). As already mentioned 
our geometrical construction must comply with duality, and this can be explicitly 
verified on an example as shown in Fig. 10 for a dominant configuration (5,15). 




Fig. 10 

(5,15) dual configurations in the ordered and disordered phase 
a ) At left, 15 bonds removed from the fully ordered phase increases the number 
of clusters by 5 and contributes to N^^(V) of Eq. (7). 
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b ) At right, 15 bonds restored from the fully disordered phase make one clus- 
ter only out of 11 sites, reducing the number of clusters by 10, and thus 
contributing to N? 15 (V) of Eq. (8), or A> 5 d 15 (V) of Eq. (9). The open dots in 
the disordered figure (at right) are the centers of the bond plaquettes. They 
coincide with sites of the dual lattice, reproducing the same figure as in the 
left. Thus the contribution to Ng 15 (V) and N£ 15 (V) of such figures are the 
same as a priori stated from duality. 



The expansion of the free energy gives similar series for the n derivative 
with respect to /3, F (n) (/3). At /3 = A 

we write the energy cumulants as 



10 

F (n) = F (n) {f3t) = {Hq) + 2z)6nQ + 2^ + ( _ } n J- C^Z m + Oiz 11 ), (13) 

for 

1 

z= vt 

where the coefficients are given in Table 4 for n varying from to 6. 

The corresponding disordered cumulants can be obtained from the ordered 
ones by use of the duality relation Eq. (5). The relation for the energies ( Ei = 
— ) is given in Eq. (3). For higher cumulants up to n = 5 we can write 



F (3, + F «3, = ^l z i [F (., + Fil)] _ _3_ [F m + F < 2)] 



q- - ( f 

<1 



(14) 

The n = (free energy) and n = 1 ( internal energy) series match the exact 
results [12] up to M = 10. The next section is devoted to a study of the series 
n > 2. 
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5. Resummation of the large q series 

We start from the result Eq. (13) which gives the cumulants of the 

ordered free energy F (f3), taken at f3 t , as series expansions to order 10 in the 
variable z = l/y/q- The coefficients of q~ m / 2 in are given in Table 4 up 
to n = 6. We want to explore the behavior in q of F^ n \ n > 2 as q is decreased 
towards q = 4. The first cumulants F ((3 t ) and F^ = —E Q are known exactly. 
For later convenience we write Fq as 

F^ = (l + z) + ^ 

where the latent heat C is given by Eq. (4) , with 9 defined through 

2 cosh 6 = y/q. 

At first glance to Table 4, the task of resumming the series p( n ) ( from now on 
we omit the index o) looks quite discouraging: not only increase all the (known) 
C^'s very fast with m, the more so n is large, but also they are all positive for 
n > 2. We will undertake this task, however, with the help of a few assumptions 
on the singularities in q of the F^'s, and after checking that the techniques used 
work well for the known case of C. 

From Eq. (4), we known that the radius of convergence of the series in q~ 1 / 2 
of C is 1/2 (q = 4). We will assume that it is so for all the F^'s. Furthermore, 
the leading singularity of C at q = 4 is given by 

^(l + ^tanh^^V 1 / 2 , (15) 



where 



7T 2 

^ = exp( — ) (16). 



A similar singularity at q — > 4 occurs in the disordered phase correlation 
length [3,4] 

Our second assumption will be that for n > 2, F^ diverges at q — > 4 like a 
power of x, up to a smooth factor. Arguments for that have been given in [8,10], 
where we proposed that i?( n )/ x 3n/ ' 2 ~ 2 is a slowly varying function of x. 
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Fig. 11 

Relative difference to the latent heat exact result of 4/6, 5/5 and 6/4 Pade's applied 
to £ _1 (dotted lines), to the logarithmic of C~ x (dashed lines) and to the regularized 
logarithmic (solid lines). The 4/6 simple Fade has a pair of pole and zero near the real 
axis giving a spurious spike near q=20. 

According to the above results and assumptions, we now experiment resum- 
mation techniques by Pade approximants on the series for which diverges as 
a; 1 / 2 , Eq. (15). Truncating it at order 10, we want to compute £ _1 from 

2/T 1 = l+z+5z 2 +7z 3 +27z 4 +41z 5 + U3z 6 +225z 7 +737z 8 + 1187z 9 +3713z 10 + - ■ ■ 

(17) 

a series which does exhibit the same qualitative aspect as F^ n \ n > 2. In the 
absence of any further information on the behavior of we would try a Pade 
resummation (P). With the knowledge we have, a better attempt is to apply the 
Pade resummation to the logarithm of £ _1 (Pin). But log£ -1 ~ 1/0 is still 
singular at q = 4; its leading singularity is a pole at 1 in the variable 

u= 2Z 2 -^=, (18) 
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Fig. 12 

Relative difference to the latent heat exact result of regularized logarithmic Fade applied 
to truncated series at 8 (dotted), 9 (dot-dashed), 10 (dashed) and 12 (solid line) terms. 

and u is expandable in z around 0. So an even better Pade technique in this case 
is to construct Pade approximants for the less singular function (1 — u) log(2£ _1 ). 
Using 2£ _1 rather than any other constant times C~ l is arbitrary, but happens to 
be very convenient as avoiding logarithms of numerical constants. Therefore we 
will construct, for any series S(z) whose lowest order term is Cfc min z fcmin and which 



is assumed to be dominantly a power of x as q 
Pade approximant as 



4, the regularized logarithmic 



(Pln R ) : S = C k 



exp 



1 — u 



Pade[(l — u) log 



Ck 



nin ^ 



(19) 



The Pade [/(z)]'s are Pm/Qn ratios of polynomials of degrees M and N, 
their Taylor expansion matches that of f(z) up to order M + N. Of course, with 
regards to the Pade techniques here applied, a series S or its inverse leads 
to the same set of approximants and the purpose of considering C~ l was just to 
illustrate its similarities with the series at hand. 
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The results for M/N = 4/6, 5/5, and 6/4 of the techniques (P), (Pin) and 
(Pln R ) applied to C are shown in Fig. 11, where the error A(£)/C = £ app ^-£ exac 
is plotted against the variable x. Typical q values are also shown. As expected, 
(Pln R ) gives the best result while (P) is worse by far (the 4/6 approximant is 
even ill behaved around q = 20). The precision reached is still less than 3% with 
(PlnR,) at q values as "small" as 5. We adopt this technique throughout the rest 
of this paper. An idea of the convergence of this resummation with the length of 
the series is given by Fig. 12 showing various results for A(C)/C using the series 
truncated at 8, 9, 10 and 12 (the 11th term is zero). The convergence is fast; the 
diagonal and near diagonal Pade's lead to comparable results. 

We now turn to the study of the F^'s with n > 2 for which the series are 
known to order 10 in z. Because their lowest term is ~ z 2 , the maximum value 
of M + N available in (Pln R ) is 8. For we find consistently that it behaves 
nearly linearly with x for x > 100 (q < 10) as well exposed by the plot of Fig. 13 
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Regularized logarithmic Fade resummations for . 

showing various estimates of F^ 2 '/x. If we assume that this quantity is indeed 
asymptotically constant, its limit is very close to 

F ( V/x ~ a = 7.610" 2 . (20) 

In any case, there is a clear cross-over around x = 100 between the large q regime 
and a very different low q — 4 behavior. Because the various Pade's lead to very 
stable results for q > 7, our result also provides accurate predictions for the specific 
heat (see below). 

Since we have a good control on F^ behavior, we next analyze F^ by the 
ratio F^/F^\ a slightly better series than that of F^ itself. With the prejudice 
that F( 3 ) ~ x 5 / 2 describes the leading behavior and with F^ ~ x, it is convenient 
to consider the positive quantity —F^/(x s ^ 2 F^) as a function of x, as shown in 
Fig. 14. The four curves correspond to the Pade's 3/5, 4/4, 5/3 and 6/2. Note 
that although F^ varies by 4 orders of magnitude in the x interval shown, the 
above ratio stays between ~ .1 and .25. If now we assume that F^ 3 ^/(F^ 2 ^ 3 / 2 ) 
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has a finite limit at large reasonable estimate is 



i?(3) 



= .17 



(21) 



x 3/2p(2) 



indicated by the dashed line in Fig. 14. 

Given the estimates a and f3 for the supposed limits of F^/x and / (x 3 / 2 F^), 
the ansatz proposed in [8,10] for the free energy of the ordered phase gives a pre- 
diction for all the higher cumulants 



We have compared direct (PfriR,) estimates of F^ n \ n = 2, • • • , 6 with the 
above parametrization where a and (3 have been fixed to their guessed values. 
The results for (-) n F^ as a function of x is shown in Fig. 15. For each n, 
four curves are drawn, one showing the ansatz (dotted lines), the three other ones 
resulting from 4/4, 5/3 (solid lines) and 6/2 (dashed lines) Pade's. The 3/5 Pade 
gives non-sense answers for n = 5 and 6, for which case the 6/2 approximant tends 
to blow up at small q. Nothing convincing can be extracted for n > 6, where the 
series become really too short. 

With this in mind, we consider the results of Fig. 15 as a manifest evidence 
for an indefinite increase with x of all the F( n )'s,n > 1 and a good indication that 
the ratios F^ /x Sn / 2 ~ 2 become smooth functions of x at x large enough, with 
relative ratios close to that of Eq. (22) (note that although F^ varies by ~ 7 
orders of magnitude when x increases from 10 to 100, its value differs by less than 
a factor 2 from Eq. (22)). 

Independently of any prejudice on the behavior of the cumulants, our analysis 
finally provides a quantitative prediction for the values of the first cumulants 
at fixed (not too small) values of q. Their knowledge may be of great help in 
understanding better the way how the thermodynamical limit is reached in the 
Potts model case, so accumulating experience on the use of finite size effects in 
other studies of phase transitions. We give in Table 5, for the ordered phase, the 
value of the specific heat at the transition point (C a = (3 2 Fo 2 ^) and of the two 
next moments F^ and F^ 4 \ together with results on C Q from existing numerical 
simulations. The uncertainties quoted contain some arbitrariness, as often when 
Pade techniques are involved. To be specific, we quote as central values the results 
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Fig. 15 

Energy cumulants (— ) n _F( n ) from 4/4, 5/3 (solid lines) and 6/2 Fade resummation and 
comparison to the ansatz Eq. (22) (dotted lines). 

from 4/4 Pade (for it is the 4/4 result of F^/F^ times the 4/4 result of 
F( 2 )). The uncertainty in each case is the mean distance to the 4/3 and 3/4 results. 
We clearly contradict values recently obtained from low temperature series [14]. 
As already noticed [8], there is strong discrepancies for q < 10 with C a values 
extracted from earlier numerical simulations [15-18]. Let us note that we quote 
two C values from [16], the value 12.7 from an analysis at the maximum of the 
specific heat and the value 18. from an analysis at the transition temperature, 
and this inconsistency was the starting point of the present work. 

In contrast our predictions agree very well with the most recently published 
values [16-21]. These data correspond to higher q value [16] or to analysis at the 
transition point in the disordered phase [20] and in the ordered phase [21]. The 
third momentum have been measured in these two last references and agree also 
very well with our predictions. 

The discrepancies are observed with data obtained from FSS analysis at the 
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maximum of the specific heat. In fact the corresponding analysis are neglecting 
energy cumulants higher than n = 2 in the expansions in 1/V as made in [15]. 
Thus the discrepancies observed can be understood from the large values of the 
higher energy cumulants as we have shown, and we can try to be more quantitative 
by introducing them in the 1/V expansions. For example the (3 value at which the 
second moment is maximum expands as 

1 Inq 1 In 2 q- 12 

PmaxF( 2 ) — Pt + 



v F (y v 2 2 F « 3 

-L^rvr[-24Ff + (Ir? q-m\nq)(F^F^ -3F™*)] 
V 6_p 

A^TTW [-H52F (1) F (2) F (3) + 1728F (2)3 
V 4 24Fi J 

+ 961n ? (4FfF! 1) F! 2) -F| 4) Fl 1)2 ) 

+ (In 4 g - 72 In 2 q)(F^ F^ - + 15Ff ')] 



+ 



where F± = F Q ±Fa and F^"^ are obtained from Fq , k = 1 to n, with the duality 
relation. At q = 10 we obtain 

3.3 2.2 1.9 10 5 1.5 10 s 

PmaxF(2) — Pt y H ^3 + • • • 

It is clear that this expansion is useless in practice as well as that of the maximum 
of the specific heat. 

Detailed comparisons between data and consequences of our results on energy 
distributions can be found in [10]. 

6. Conclusions 

We have explored the properties of the 2D-Potts model free energy F{j3) in the 
region q > 4 where the model has a temperature driven first order transition. This 
was achieved by a series expansion of F(f3), close to the transition temperature, in 
powers of 1/y/q, performed to order 10 from the Fortuin-Kasteleyn representation 
of the partition function. At each order m in q~ 1 / 2 , we compute the number of 
lattice configurations of / bonds enclosing k clusters of sites such that / — 2k = m. 

The results obtained were translated into similar series in q~ x l 2 for the deriva- 
tives F^ of the free energy taken at the transition temperature, that is for the 
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energy cumulants of the model. These series truncated at order 10 provide strong 
evidence for highly divergent behaviors at a low q value, the more so n increases. 
We assumed that this value is q = 4, and conjectured that these singular behav- 
iors are dominantly embodied in powers of the variable x(q) of Eq. (16), known 
to accurately describe the latent heat and the largest correlation length in a wide 
region of q > 4. Then, Pade techniques adapted to such behaviors were applied to 
the series, leading to quantitative predictions for n < 4 , q > 6, confirming 
severe divergences of all F^'s and in semiquantitative agreement, for n up to 6 
at least, with an ansatz [8-10] which prescribes the divergence rates. 

First these results enlarge our analytical knowledge of the 2D-Potts model, so 
providing useful additional tests on methods for analyzing numerical data on finite 
lattices. Next they illustrate the possible occurrence in first order transitions of 
properties widely influenced by a nearby continuous transition point. This may be 
the case as well for field driven phase transitions just below the critical temperature 
when the correlation is still large [22,23]. 

One may think of other applications of our analytical approach to the free 
energy of the Potts model. An interesting ones would be a study of the analytic 
structure of F (f3) ( or Fd(f3) ) at fixed q > 4 around (3 = f3 t , where an essential 
singularity is expected [24,25]. 

We thank A. Billoire and P. Moussa for useful discussions and N. Elstner for 
a critical reading of our manuscript. 
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Table 1 

Number N of unequivalent (upon lattice translations) dominant configurations. 
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Contribution 


Cluster set 




+2 C(7 5) 


(2 7) 


i 


-C2y - 16) ecu i) 


(2 7) * (1 4) 


2 


-12 C(10, 2) 


((2 7) + 3-3 10) 


3 


-4 C(12, 1) 


r r 2 7) + 4-4 12) 

VV ' / 1 > * / 


4 


-4 


CC2 71 + 5-4 12) 

VV > / ' ' ' / 


5 


-28 


CC2 71 + 5-4 13) 

VV 5 1 J ' °1 1 "J 


Q 


+8 


(((2 71 + 3-3 10) +2-4 12) 

VVV > / 1 ' / 1 > ' / 


7 


+12 CflO 3) 


(2 9-3 10) 


g 


-8 CC12 1) 


(72 9- 3 10) + 2-4 12) 


9 


-80 


((2 9- 3 10) + 3-4 13) 


10 


+6 Cfl2 2) 


f2 10-4 12) 


11 


+48 C(13 1) 

| J- \ ' I _L '.J % -3-1 


(2 11-4 13) 


12 


+8 C(15, 1) 


(2, 11; 5, 15) 


13 


+80 


(2 12- 5 15) 


14 


+12 


(2 12-6 17) 


15 


+1/2(V - 5) C(8 4) 

-1/2(V 2 - 15V + 62) 


d 4) * ci 4) 


16 


d 4) * d 4) * ci 4) 


17 


-12/2 CflO 2) 


((1 4) * fl 4) + 2- 3 10) 

VV > / V > / ^ 7 5 / 


18 


-4/2 


((1 4) * fl 4) +4-4 12) 


19 


-52/2 


((1 4) * d 4) +4-4 13) 


20 


+8/2 


CCCl 4)*fl 4) + 2- 3 10) + 2- 4 12) 


21 


-2/2(4V-32) C(ll, 1) 


(1,4)* ((1,4) + 3; 2, 7) 


22 


+ (2V - 16) C(ll,2) 


(1,4)* (1,6; 2, 7) 


23 


-28 


((1,4)* (1,6; 2, 7) + 2; 4, 13) 


24 


+ (6V -62) 


(1,4)* (1,8; 3, 10) 


25 


+ (V-12) 


(1,4)* (1,8; 4, 12) 


26 


+l/2(4V-46) 


(1,6; 2, 7)* (1,6; 2, 7) 


27 



Table 2 

Detailed contributions to ^2,12/V ^ column 1 (with C(i,j) as 
given in Eq. (10) ) of the different cluster sets in column 2. 
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48 
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1 00263 


23 






14 


-6054 


22 




6 


1182 


-379348 


21 






-244 


71360/3 


20 




22 
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3321646/5 


19 






-208 


-78920 


18 




6 


6704 


-625246 


17 




2 


-138 


106586 


16 




-63 


-65917/4 


351774 


15 




8 


2156 


-226216/3 


14 




-201 


15532 


-119493 


13 




6 


-3000 


29954 


12 


1 


1555/3 


-7365 


21723 


11 




-72 


1762 


-6066 


10 


6 


-406 


1686 


-6431/5 


9 




88 


-468 


3422/9 


8 


-33/2 


131 


-459/4 
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2 


-38 


254/7 




6 
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-37/3 
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-4 


22/5 
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1 -3/2 
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Table 3 

Non zero coefficients A(l,m) contributing to F D in Eq. (12) 
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m\n 





1 


2 


3 


4 


5 


6 


2 


o 


4 


16 


64 


256 


1024 


4096 


3 


-4/3 


-2 


34 


430 


3778 


29518 


218914 


4 


1 


2 


114 


2654 


41778 


556382 


6813714 


5 


-8/5 


-6 


254 


12186 


322670 


6773994 


126069374 


6 


2 


4 


882 


57018 


2210982 


67122114 


1774583142 


7 


-12/7 


-16 


1944 


224732 


12819264 


546094604 


19774354944 


8 


5/2 


6 


6128 


888024 


68657204 


3918393456 


187361651588 


9 


-4/9 


-38 


13550 


3164682 


333583598 


25037212842 


1545876302510 


10 








39698 


11243178 


1532324246 


146961943266 


11451708807878 



Table 4 

(n) 

Coefficients of the expansion of F in Eq. (13). 



Q 


Co 


_ F ( 3 ) 

■T o 


-To 


^exp 


_^(3)exp 


Ref. 


30 


3.41294(5) 


16.74(2) 


7.00(5) 10 2 








20 


5.3612(4) 


56.9(4) 


5.0(1) 10 3 


5.2(2) 
5.38(4) 
5.351(15) 


55.8(9) 
57.0(13) 


[19] 
[20] 
[21] 


15 


7.999(3) 


179.(4) 


3.1(2) 10 4 


8.04(4) 
8.016(21) 


175(5) 
180.5(31) 


[20] 
[21] 


10 


17.98(2) 


1.9(2)10 3 


1.3(2) 10 6 


10.6(11) 
12.7(3) 
~18. 
18.0(2) 

17.95(13) 


2066(81) 
1979(87) 


[15] 
[16] 
[16] 
[20] 
[21] 


8 


36.9(2) 


1.4(4)10 4 


2.7(8) 10 7 


22.8(30) 




[15] 


7 


69.6(5) 


7.(3)10 4 


3.4(13) 10 8 


47.5(25) 
50.(10) 
44.(22) 




[17] 
[16] 
[18] 



Table 5 

Results for the first energy cumulants at some q values 
and comparison to numerical data. 
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